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We derive the connection between the Cooperon problem in weak localization theory and the random 
matrix description of type-II superconductors. As magnetic field and disorder increase, an extreme 
type-II superconductor crosses over from a state in which the low energy quasiparticles are primarily 
localized and the density of states is determined by the electronic structure of individual vortices, 
to a 'chaotic' state, in which quasiparticles are primarily extended and the density of states is 
determined by the random matrix description. 
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I. INTRODUCTION 

In this paper we explore the hypothesis that electrons 
in a superconductor with explicit time-reversal breaking 
pair via random matrix elements. This hypothesis com- 
bines ideas and techniques from three fields of study usu- 
ally considered unrelated: random matrices, type-II su- 
perconductors, and the Cooperon problem of weak local- 
ization. We will show that incorporating these ideas into 
a random matrix model allows us to calculate the density 
of states of a type-II superconductor non-perturbatively 
and leads to new possibilities for investigating chaotic 
behavior in complex quantum systems. 

This section contains brief introductions, with refer- 
ences to extensive reviews, of the main ingredients of 
the theory. In section II we describe the random ma- 
trix model, a generalization of the Anderson description 
of disordered superconductors to superconductors in a 
magnetic field. In section III we derive the connection 
between this description of type-II superconductors and 
the Cooperon. In section IV we show that this connec- 
tion yields the correct answer for the semiclassical upper 
critical field, Hc2- In section V we discuss the weak field 
limit and the relation between the random matrix model 
and the Abrikosov-Gor'kov theory. Finally, in section 
VI we summarize the principal results and describe open 
questions and further possible areas of research. 



A. Random Matrices 

Random matrices were introduced into theoretical 
physics by Wigner in the 1950's as a tool for under- 
standing the distribution of energy level spacings and 
widths observed in nuclei [||. The distribution of level 
spacings in nuclei was compared with the distribution of 
level spacings of a large matrix filled with random (un- 
correlated) matrix elements; quantitatively good agree- 
ment was eventually found . The idea of comparing the 
level spacings of complex quantum systems with the level 
spacings of random matrices was eventually extended 



to atoms H], small metallic particles [Q, and classically 
chaotic systems [§. In all of these applications, in con- 
densed matter, atomic, and nuclear physics, the focus has 
been on understanding the distribution of the spacings of 
energy levels, that is, the probability P{E,x)dx that the 
neighbor of a level at energy E has an energy in the range 
E + X to E + X + dx. In contrast, we will use random 
matrices to understand the density of energy levels, that 
is, the number of energy levels p{E)dE in the range E to 
E + dE. 

Random matrices have also been studied in the con- 
text of field theory, spurred by 't Hooft and Brezin et 
al.'s discovery of the connection with the large-N limit 
of certain SU(N)-invariant field theories A connec- 
tion between the Feynmann diagram expansion for ran- 
dom matrix theories and the sum over random surfaces 
in two-dimensional quantum gravity was also established 
0. More recently, random matrices have resurfaced as 
a possible unification of various string theories and 
as a concrete example of a more general theory of "non- 
commuting random variables" |^ . 

Extensive recent reviews can be found in Refs. [ pO| and 



B. Type-II Superconductors 

The BCS theory [|l2| describes the electronic properties 
of ideal, uniform superconductors; it applies to the interi- 
ors of type-I superconductors, which expel magnetic field, 
or to type-II superconductors below the field of first pen- 
etration, Hci- The presence of a magnetic field, which 
occurs for type-II superconductors in the mixed state, 
Hci < H < Hc2, makes the calculation of electronic 
properties significantly more complicated. The magnetic 
field destroys both the spatial uniformity and the time- 
reversal symmetry implicit in the BCS description. 

Magnetic field penetrates a superconductor in the form 
of tubes of magnetic flux (vortices) enclosing circulat- 
ing quasiparticle currents. These currents lead to elec- 
tronic states which are bound to individual vortices. Such 
states are not only spatially inhomogeneous, but may 
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overlap strongly with neighboring vortices and form en- 
ergy bands, the details of which will be sensitive to the 
presence of disorder or the deviations from an ideal vor- 
tex lattice that are found in most materials. 

It is possible to attack the problem of the electronic 
properties of type-II superconductors numerically. For a 
fixed type of vortex lattice, a large unit cell can be con- 
sidered which contains one or more vortices. Either the 
full Bogoliubov-de Gennes equations or the approx- 
imate Eilenberger equations (derived by assuming 
slow variations on the scale of the Fermi wavelength) may 
be solved to obtain the quasiparticle eigenstates [ p^|p^ . 

The numerical approach to the electronic properties of 
type-II superconductors has several drawbacks: (1) It is 
difhcult to get useful results numerically, both because of 
the large number of electrons per unit cell for realistic ma- 
terials and because of the widely different energy scales in 
the problem. (The cyclotron frequency is much less than 
the BCS gap, which is in turn much less than the Debye 
energy.) (2) In a magnetic field, unlike conventional su- 
perconductors at zero field, the electronic properties are 
sensitive to disorder. The highly idealized impurity-free 
perfect vortex lattice is unrealistic for most materials. 
(3) The numerical approach calculates every quasipar- 
ticle wavefunction in the superconductor. This is often 
an overkill: for many applications we need much less in- 
formation, for example, just the single-particle density 
of states. An approach which is simpler and more gen- 
eral than numerically solving the Bogoliubov-de Gennes 
equations is desirable. 

Analytically manageable extensions of the BCS theory 
to disordered superconductors were given both by An- 
derson [|r^ and by Abrikosov and Gor'kov [Q. Ander- 
son generalized the BCS description to show that non- 
magnetic impurities have little effect on conventional su- 
perconductors. Abrikosov and Gor'kov introduced a di- 
agrammatic technique to include the effects of a random 
ensemble of impurities, and showed that magnetic impu- 
rities can suppress superconductivity. 

Neither the Anderson nor the Abrikosov-Gor'kov de- 
scription, however, applies directly to the mixed state 
of type-II superconductors. The Abrikosov-Gor'kov de- 
scription requires spatial homogeneity, which is violated 
in the mixed state, and, in order to apply to a type-II 
superconductor, must assume that random magnetic im- 
purities are equivalent to an applied magnetic field [p^ . 
The Anderson description is more general, in the sense 
that it is not tied to the assumption of averaging over a 
random ensemble of impurities, however it requires both 
time-reversal symmetry and spatial homogeneity. The 
method described below can be considered a generaliza- 
tion of the Anderson description which has neither of 
these constraints. 



C. The Cooperon 

The problem of the localization of a quantum particle 
by a random potential was originally posed by Ander- 
son in 1958 ||2^, in the context of understanding metal- 
insulator transitions, and has led to a large body of work 
on the subject of the quantum effects of disorder in met- 
als. The field has been summarized in a number of recent 
reviews pl|-p^. 

One of the theoretical approaches to understanding the 
quantum effects of disorder is impurity-averaged pertur- 
bation theory . For non-interacting electrons scat- 
tered by rigid impurities. Green's functions may be ex- 
panded perturbatively in the impurity potential. Averag- 
ing over a random ensemble of impurity potentials yields 
classes of diagrams which can be organized by the powers 
of kp£ they contain, where kp is the Fermi momentum 
and i is the mean free path. Diagrams within a given 
class may be summed, resulting in a series expansion in 
kp£ for the desired Green's function. 

Impurity-averaged perturbation theory can be used to 
calculate the conductivity of an electron gas in the pres- 
ence of random disorder. The conductivity is related 
to the current-current response function by the Kubo- 
Greenwood formula, and the impurity diagrams for the 
current-current response function can be summed. It was 
recognized early on that a certain class of diagrams, the 
maximally crossed ones, yield a significant quantum cor- 
rection to the classical conductivity This estimate 
for the quantum correction was used as support for the 
scaling theory of localization |^,^ . 

The maximally crossed diagrams for the current-cur- 
rent response function (a particle-hole response function) 
can be related by time-reversal to ladder diagrams for a 
response function in the particle-particle channel. The 
impurity-averaged particle-particle response function was 
called the 'Cooperon' by Altshuler et al. | |2^ , p0| . It is this 
response function, the Cooperon, originally used to study 
the conductivity of a disordered electron gas, which ap- 
pears in the random matrix description of type-II super- 
conductors described below. 



II. RANDOM MATRIX MODEL 

The starting point for an electronic description of type- 
II superconductors is the Hamiltonian |^ 

n = Jdrcl^noir)c^a~^Von j dvcl„cl^,c^„,c^„. (1) 

Here 51 is the system volume, Vo > represents an at- 
tractive short-ranged interaction, the creation operators 
satisfy 
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and the bare Hamiltonian Tip is 



V-^A(r))' + C/(r)-i?j. , (2) 



with U{r) the potential due to a static distribution of im- 
purities. The creation operators in the interaction term 
in Eq. are understood to be superpositions of opera- 
tors which create states only near the Fermi surface, that 
is, Cr = Jdp exp{«p • r)cp0p, where the theta-function 
limits the energy to be within a Debye frequency of the 



Fermi energy: 9p = 9{uJi 



Ep\). In the weak- 



coupling limit, where the range of the attractive poten- 
tial is much smaller than the extent of a Cooper pair 
{vp/ujD ^^v/Ao), the interaction term may be consid- 
ered local p2| . 

The variational approach to the Hamiltonian in Eq. 
(0) is defined by 



W ^ dr ^ 



Hoiv) $(r) 
a>*(r) -H5(r) 



(3) 



where 



The order parameter $(r) is determined self-consistently 
from the equation 



$(r) = -yo^^(CrTCr^) 



(4) 



For a pure system {U — 0) with no magnetic field, plane 
waves diagonalize the bare Hamiltonian Ti.Q. If the or- 
der parameter is assumed to be a constant, $(r) — Aq, 
then Eq. separates into 2x2 matrices. This yields 
the standard BCS results: the spectrum is given by 



where Ek 



Ep, and Eq. 



becomes the BCS gap equation pqVq J de]^/2Ek — 1, 
with po the density of states at the Fermi level. 

In the case of a system with impurities (J7 7^ 0), but 
with no magnetic field, a similar result holds, as shown 
by Anderson . If we write the eigenstates of the bare 
Hamiltonian as 



and rotate to a basis defined by these eigenstates, 

dv V'°(r) 4, , 
then the off-diagonal, pairing term in Ti' becomes 
j dv $(r) = A ^ J^^f, 

a/3 



1/3 7 



(5) 



(6) 



(7) 



where 



dv i,l{j)^%{v) 



(8) 



Note that we explicitly assume that the order parameter 
is spatially homogeneous: ^(r) = A. The time- reversal 
symmetry of the system (no magnetic field, no magnetic 
impurities) implies that Tig = and therefore a basis 
may be chosen in which i/jJ^ = 'i/'a*- The overlap integral 
in Eq. (ph then yields 



(9) 



The Hamiltonian H! again factorizes into 2x2 matrices, 

the spectrum is again given by E^ — \J e1^ + h? and solv- 
ing the self-consistency equation yields the same result as 
in the impurity free case (assuming the density of states 
at the Fermi level is unchanged): A — Aq. 

The random matrix description we will use for type-II 
superconductors is a generalization of the Ander- 
son description in the sense that we start by using ex- 
act eigenstates of the bare Hamiltonian TLq as a basis. 
This includes whatever vortex or impurity distribution is 
present. The eigenstates are 



and we rotate to the basis 



dv V'a(r) Cr 



(10) 



(11) 



There are two important differences from the zero mag- 
netic field case. First, the order parameter will no 
longer be spatially homogeneous; the phase winds rapidly 
throughout the superconductor (by 27r near the center of 
each vortex) and the magnitude vanishes at the center of 
each vortex. In order to extract a typical measure of the 
magnitude of the superconductivity, we define 



$(r) EE 0x(r) 



(12) 



where x(r) is normalized so that / |x(i")P dv — Q.. Hence 
(fp' is the spatial average of the magnitude of the order 
parameter. Second, we no longer have time-reversal sym- 
metry, because of the applied magnetic field. We still 
define the pairing matrix as 



A^p = ydr x(r)Va(r)^0(r) , (13) 
but we no longer have Aafi = Sap- Specifically, writing 



Aaf3 = gap h^^'^iea - ep) , 



(14) 



then in zero magnetic field we have h{sa — £p) = 5ap- 
In the next section we will show that for a randomly 
disordered electron gas in a magnetic field, the i5-function 
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distribution for h{s) broadens by an amount proportional 
to the magnetic field. 

The motivation behind the random matrix model is 
that ga/3 is a material-specific, rapidly varying complex 
function, the details of which should not affect the den- 
sity of states averaged over a random ensemble of dis- 
ordered superconductors. The average density of states 
is sensitive to the overall structure of the pairing matrix 
elements ~ the broadening of the zero-field ^-function dis- 
tribution of h(e) - but not to the rapidly varying 0(1) 
complex numbers which multiply it. 

The basic assumption of the random matrix model is 
therefore that average properties of an ensemble of super- 
conductors can be modeled by choosing complex numbers 
gafs from a random, uncorrelated distribution: 



9al3 9a' p' 



(15) 



where the brackets denote an average over the distribu- 
tion: [X]av = (l/iV^)EM=i^^''^- The particular choice 
of symmetry for ga^5^ whether real symmetric (GOE), 
Hermitian (GUE), or complex symmetric, will affect the 
eigenvalue spacing distribution but will not affect the 
density of states in the large matrix size limit. 

Using the above definitions for <i>(r) and A, and rotat- 
ing to the basis of bare eigenstates defined by Eq. (11), 
the Hamiltonian in Eq. (H) becomes 



Eo (t)A 



(16) 



where Eq — diag(. ..£„...) is the diagonal matrix of the 

eigenvalues of Ho, and ^I^t = [. . . c^.^ cpi . . .]. The 

assumption of uncorrelated random numbers for gap al- 
lows us to determine exactly, in the large matrix limit, 
the Green's functions of H' [BSl . 



III. THE PAIRING MATRIX AND THE 
COOPERON 

The random matrix assumption for gap lets us write 
down integral equations for the Green's functions of the 
Hamiltonian Ti' for any value of 0, the overall strength 
of the pairing, and h{uj), the average pairing amplitude 
defined in Eq. (p^. This will only be useful in comparing 
with experiment if we have some prediction for, or theo- 
retical understanding of, the pairing amplitude h(uj). For 
a diffusive system, and a disordered vortex lattice, there 
is a simple prediction. 

We can rewrite Eq. ( p^ ) to express h{Lu) in terms of the 
average squared matrix element as a function of energy 
difference: 



(17) 



where the brackets denote averaging over the random 
ensemble. The density of levels of the normal metal, 
Pq = p{Ep), is assumed to be independent of energy 
in the narrow range about the Fermi energy relevant to 
superconductivity. The total number of levels that are 
being paired in Ti' , which determines the size N of the 
matrix, is the level density times the maximum pairing 
energy. This cutoff is usually taken to be the Debye en- 
ergy, so = pQiOu- The normalization factor l/Npo in 
Eq. (|l^) follows from Eq. (14) and the continuous limit 
for the energy levels: J2a d£ap{£a)- 

Note also that we will always be interested in the limit 
where N oo, because is much larger than the level 
spacing (as well as the pairing strength (j) and the bare 
BCS gap Aq). In order for the Hamiltonian H' in Eq. 
( |l6| ) to be interesting in that limit, the matrix elements 
Aai3 must be 0{1/Vn). (This is easy to see with per- 
turbation theory: the shift in energy of some level is 
Ssa = </'^S/3 |-^a/3p/(£Q + ^ p) ■ The TV terms in the sum 
must be canceled by a 1/N factor from |^Q,^p.) 

Inserting Eq. (|l|) for Aap into Eq. ^ yields 



X x(r)x*(r') 5{u:-ea+ef)) 



(18) 



We can insert a factor of 1 = / dE S{E — Sa) and write 
this as 



Utj) = Idrdr'dE 
Npa 



6Go{r,r',E) 



X 5Go(r,r',^ + L^) X(r)x*(r') , (19) 

J av 

where the bare single particle Green's functions are 

SGq = Gq — Gq 

MrWair') 



G^{r,r',E)^J2 



To make further progress, we will need to make a key 
assumption: that there is, on average, no correlation be- 
tween the bare single particle Green's functions G^ and 
the order parameter product x(r)x*(r') appearing above: 



G^G^XX* 



*-^o 



XX 



(20) 



The wavefunctions which enter Go are the eigenstates of 
the bare Hamiltonian Hq, not the eigenstates of the full 
superconducting Hamiltonian. Hence the assumption is 
that the product of the bare Green's functions is uncorre- 
lated, on average, with the order parameter, that is, the 
positions of the vortices. Below, we will only be inter- 
ested in the case in which the magnetic field is nearly uni- 
form within the superconductor, H 3> Hd- In this limit, 
since the magnetic field in Ti.Q is uniform, the positions 
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of the vortices affect Go only through a gauge-dependent 
phase. The quantity we are interested in, h{u!), is gauge- 
invariant. Therefore, with the gauge- invariance of the 
product of the two terms properly maintained, and for 
H ^ Hci, there are no additional correlations between 
the bare Green's functions and the order parameter. 

With the average over the order parameter product 
written as 



X(r)x*(r') = g(r,r') , 

J av 

the expression for h{uj) then becomes 



(21) 



h{uj) 



(2^ 

Npo 



Re / dr dr' 



dE G3"(r, r',E) 



X Goir,r',E + u;) g{r,r') , (22) 

- av 

where we have used the identity 



(G+{r,r',E) {r,r' , E + u;))* ^ 



Go{r',r,E) G+ {r' ,r, E + cj) 



(23) 



and we omit the terms involving G^G^ and G^G^ be- 
cause they do not contribute. 

As mentioned in the introduction, in weak localiza- 
tion theory quantum corrections to classical transport 
are calculated by studying averaged values of response 
functions of non-interacting electrons in a random disor- 
der potential. The quantity in brackets in Eq. (|2^) is one 
such disorder-averaged response function: the Cooperon 

HQ- 

The order of r and r' in Eq. ( p2| ) is important. If 
the single particle Green's functions appeared in the 
form Go(r, r', E) Go(r', r,E + uj), the quantity in brack- 
ets would be just the real-space Fourier transform of the 
usual density-density response function. It is because 
Aap is a matrix element for electrons in a Cooper pair 
(the overlap between Tpa and ^pp rather than t/;* and ipfs), 
that the order of r and r' is reversed, and the bracketed 
quantity in Eq. ( p2| ) is the Cooperon rather than the 
diffuson. 

The Cooperon response function G(r, r',w) is defined 

as 



G(r, r',uj) 



dEG+{v,v',E) G^{r,r',E + iu) 



(24) 



This response function, in the limit in which the mean 
free path is much smaller than the magnetic length but 
much larger than the Fermi wavelength (Af ^ ^ ^ = 
y^hc/eH), obeys a diffusive equation of motion |2|] 

2 



-iuj + D [iV + I^A 

he 



G(r,r',w) = 5(r- r') . (25) 



The diffusion constant D = Vp£/d, where d is the space 
dimensionality. This equation defines the Green's func- 
tion for the Schroedinger equation of a single particle of 
charge 2e in a magnetic field. The solution is therefore 



G(r,r',c.) = 



27r -iuj+{AD/ijj){n + 1/2) + Dk^^ 



(26) 



where ifnkix,y) are normalized 2D Landau level wave- 
functions and we assume a uniform magnetic field H 
(vahd for H > Hci). 

This solution is not sufficient to evaluate h{uj); we also 
need to know the order parameter correlation, g(r,r'), 
defined by Eqs. (|l^) and (pi]). The phase of the order 
parameter winds by 27r around each vortex. In a su- 
perconductor with a perfectly periodic vortex lattice we 
expect $(r)$*(r') to refiect the long-range ordered struc- 
ture of the vortex lattice. Averaged over an ensemble 
of disordered materials, however, the phase will disorder 
on the length scale set by the inter- vortex spacing 
Therefore we will assume that the averaged correlation 
function has the form 



5(r,r') = e 



o(r,r') 



(27) 



The gauge-dependent phase factor is required to keep the 
product G(r, r' ,uj)g{r, r') gauge-invariant. In the Landau 
gauge A = Hxy, we have 6'gaugo(r, r') = (a; -I- x'){y - 

The form of \g{r, r')| is an additional assumption that 
is necessary for evaluating h(Lu). It is not essential that 
this form be Gaussian. However, the Gaussian form is 
natural for a random system, and also has the property 
that the result for h{u!) is particularly simple. 

We now substitute the expressions for C{r,r',uj) and 
g(r,r'), Eqs. (|6|) and (|3), into Eq. (H). The integral 
over z ^ z' sets kz — in the denominator. The Gaus- 
sian form for g(r,r') can be written as \g\ — ^/nhQ{{x — 
x')/£f^)hQ{{y — y')/^_f/)j where hn(u) are the normalized 
harmonic oscillator eigenstates (Hermite functions). The 
2D Landau level wavefunctions in Landau gauge have 
the form (pnk{x,y) = A exp{ikx) hn{kifj/^/2 + ^/2y/i^). 
Rewriting the product in Eq. ( p6| ) using the identity 
K,ia)hn{b) = CLh2n-jiia + b)/^/2)h,{{a - b)/V2) 
, we see that the spatial integrals in Eq. (|2^) and the 
orthogonality of the Hermite functions set n — 0. The 
final result is therefore that only the vanishing n and kz 
terms survive in the sums in Eq. (|2^). This yields for 
hiu): 



h{uj) 



W 



2npo uo^ + W^ ' 



where the width W is 



W = 2hD/£\j = 2eDH/c 



(28) 



(29) 
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and we have normahzed to po J dco h{uj) = 1 . 

The Lorentzian form for h{uj), Eq. (^), is what might 
be expected on general grounds in a diffusive system 
where the time-evolution of pair correlations obeys an 
exponential decay law. This point was discussed by de 
Gennes in the context of the "ergodic evolution of the 
time-reversal operator" |Q. de Gennes considered the 
problem of evaluating the transition temperature Tc{H) 
for a general superconducting system in a magnetic field 
described by the variational Hamiltonian in Eq. (Q). He 
used a perturbative expansion in the order parameter 
$(r), which is valid near the phase transition, and found 
that the average over single-particle wavefunctions which 
appears in Eq. is the kernel operator acting on <l>(r) 
in the gap equation, de Gennes related that average over 
single-particle wavefunctions to the time evolution of the 
time-reversal operator, K{t), for one electron moving in 
the potential described by the Hamiltonian T-L^ (Eq. (||)). 
An "ergodic" system was defined as one in which the 
motion of the electron is diffusive, in the sense that the 
time-reversal operator decays exponentially with time: 



lim {K\Q)K{t)) = e-*/^ . 



(30) 



For a system defined as ergodic in this sense, ft.(w), which 
is the power spectrum of the operator K{t), has the 
Lorentzian form given in Eq. (^), with W = t^^ . 

The interesting result here is that the average over 
the single-particle wavefunctions in Eq. (|8|) appears not 
only in the equation for Tc{H), but also in the random 
matrix description of the pairing Hamiltonian. The func- 
tional form of h{u}) and the random matrix technique al- 
low us to go beyond perturbation theory and therefore 
beyond calculating Tc{H). At any field for which the 
random matrix element hypothesis is valid, we can ob- 
tain the full Green's functions non-perturbatively. 

We note also that the connection with "ergodic super- 
conductors" suggests that the random matrix description 
may apply, in addition to the bulk disordered systems 
considered here, to other systems discussed by de Gennes 
et al. as being ergodic, for example, thin films and su- 
perconducting nanoparticles p5|,p6[. 



IV. UPPER CRITICAL FIELD 



Using Eqs. ( Pq ) and (|2^) we may derive a relation be- 
tween the upper critical field Hc2 and the diffusion con- 
stant D. This will allow us to test whether the formalism 
described in the previous section yields a result consistent 
with the prediction of the scmiclassical theory of type-II 
superconductors. That theory is based on the pertur- 
bative solution of the Gor'kov equations near Hc2 using 
semiclassical Green's functions to include the effects of 
the magnetic field. 



The difference in energy between the normal metal 
ground state of the Hamiltonian in Eq. (^, given by 
(f> = 0, and the superconducting ground state, is, to sec- 
ond order in (t>: 



AE 



--E 



af3 



0(04) . (31) 



The I^Iq/jP term is second order perturbation theory for 
the Hamiltonian Eq. (||), where 9{x) is the step func- 
tion (for non-zero temperatures, 9{eaSf}) is replaced by 
[tanh(£a/2r)-|-tanh(£/3/2r)]/2) and the upper cutoflFs for 
the energies are \ea\ < wd, l^ffl < ^d- The 1/Vb term is 
the usual Hartree-Fock self-consistency contribution and 
follows from (J dr <i>(r)cr|Crx) and $(r) — — ilVo(£rTCri)- 
We must ignore small non-perturbative effects (the 
breakdown of second-order perturbation theory) , in order 
to reproduce the semiclassical limit. 

The upper critical field is the field Hc2 at which the 
quantity in brackets in Eq. ( |3T| ) vanishes. At lower fields, 
the coefficient of (jy^ is negative and > minimizes the 
energy; at higher fields, the coefficient of 0^ is positive 
and = minimizes the energy. 

Averaging over the ensemble amounts to replacing 
|-4a/3p in Eq. (|3l|) by its average value h{sa—£i3)- Taking 



the continuum limit for the energies, inserting Eq. (|2J 
for h(h'), and using the definition of the BCS coupling 
constant, 1/g = 1/pqVq, yields that the 0{(j>'^) term in 
Eq. ( pl| ) vanishes when 



W 



£a + e/3 (Sa " Sfif -f 



(32) 



Performing the integrals yields 



1 , /2cjD 
- = log 



1 



In the weak coupling limit Wjio-D 
zero-field BCS result, 



(33) 

0. Subtracting the 
(34) 



where Aq is the zero-field BCS gap, we obtain the equa- 
tion which determines the upper critical field: 



Ao 



Inserting Eq. (C 



for W{H) yields 
27r2 Ho ' 



(35) 



(36) 



where = hvp/nAo is the zero- field BCS coherence 
length and D = vpi/i in three dimensions. 

This result may be compared with the result derived 
with the semiclassical theory of dirty type-II supercon- 
ductors. The semiclassical approach starts with plane 
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wave electronic eigenstates, assumes a weak magnetic 
field, and includes the eS'ects of impurities through the 
Abrikosov-Gor'kov impurity averaged perturbation the- 
ory technique. The resulting Gor'kov equations can be 
expanded perturbatively in the order parameter (near 
Hc2) and yield the equation which determines Hc2 1 38 - 4^] 



In 



2<^oT 



(37) 



where ipi^) is the digamma function, T^q is the zero- field 
critical temperature, and denotes the dirty limit, 
e « ^0- As r ^ this gives = (3/27r2)($o/^eo), 
in agreement with Eq. (^). 

The agreement with the semiclassical theory shows 
that the random matrix description combined with the 
results from the Cooperon problem reproduces the per- 
turbative result of the more conventional approach to 
understanding dirty type-II superconductors. 



V. CONNECTION WITH ABRIKOSOV-GOR'KOV 
THEORY, AND THE H ^0 LIMIT 

As pointed out in Ref. in the weak field limit, for 
which the T-breaking scale W ^ Aq, the random matrix 
model described in the previous sections does not yield 
the density of states the zero-field BCS theory. What 
happens in this limit is that the central assumption of the 
model, that the Hamiltonian of the superconductor can 
be well described by a large matrix with random pairing 
elements, breaks down. For small W/Aq there exists a 
special basis, discussed below, in which the Hamiltonian 
has a simpler form. The existence of a special basis means 
that correlations have developed which are not obvious 
in the original basis. (As a simple analogy, consider an 
electron in a cylindrically symmetric potential. Random 
matrix elements in a plane wave basis would be a poor 
assumption.) 

To understand what happens to the random matrix 
model of the previous sections in the weak field limit, we 
first write down two additional possible random matrix 
models for superconductors with T-breaking. The first 
of these turns out to be a random matrix formulation of 
the Abrikosov-Gor'kov theory. The second of these we 
will argue describes the weak field limit, H <C i?c2, of a 
superconductor in a magnetic field. 

The theory described in the previous sections we de- 
note by Hj -. 



Hi = 



with A satisfying 



Eq 4>A 



Sik Sji h{ei - Ej) . 



(38) 



As before, Eq = diag(. . . Eq. . . .) is a diagonal matrix of 
uniformly distributed eigenvalues. 

A different possible random matrix model occurs if we 
start with BCS and add a random matrix to the diago- 
nal components, in a way which breaks the time-reversal 
symmetry: 



H 



II 



Eq 


A 


+ a 


■ M 


" 


A 


-Eq_ 





M 



(39) 



with M satisfying 



kl\ 



Sik s 



k Oji 



Eq is a diagonal matrix as in the previous case, A is 
proportional to the identity matrix: A = Ao(5ij, a is 
a parameter which measures the amount of T-breaking, 
and M is a random Hermitian matrix with no weighting 
factor h{ei—ej) in it. The T-breaking in this model is rep- 
resented by the same sign for M in the upper and lower 
diagonal quadrants. Since these quadrants correspond to 
the energies for up-spin electrons and down-spin holes, 
this indicates a potential which has opposite sign for up 
and down spins. In other words, this model appears to 
correspond to random magnetic impurities. In fact, the 
equations for the Green's functions of this random matrix 
model are the same as the Abrikosov-Gor'kov equations 
for a superconductor in the presence of random magnetic 
impurities. This random matrix model is a reformulation 
of the AG theory H. 

The final random matrix model we consider is 



- 



Eq 


A 




A 


-Eq 


+ ax 



012 





M(2) 
M(2) 







(40) 



with the matrices M^°-'> satisfying 



This differs from Hjj in that a random matrix appears 
in both the diagonal and off-diagonal quadrants. 

How do wc best describe a superconductor in the weak 
field limit. El <^ i?c2? The order parameter in a super- 
conductor may be written as 



Hr) = /(r) e 



z9{r) 



(41) 



where /(r) is real and d{r) is a position dependent phase 
that winds by 2t: around the center of each vortex. For 
a given material, with some order parameter phase re- 
alization 0{r), we may transform the standard creation 
operators to a new basis defined by 



r ^ r P*S(r)/2 



In this basis the pairing interaction takes the form 



(42) 
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$(r) cl 



/(r)A 



(43) 



The effective order parameter in the new basis is purely 
real. Note that although this has sometimes been called 
a gauge transformation, it is not: the magnetic field 
changes. The vector potential is A(r) = A(r) + V6'/2e so 
/ A(r) • dl ^ / A(r) • d\ around a loop enclosing a vortex. 

Physically, we are transforming to a basis in which the 
magnetic field is the previous field distribution plus one 
flux quantum threaded through the center of each vor- 
tex in a direction opposite to that of the applied field. 
This means that far away from an isolated vortex, the 
net effective flux seen by an electron is zero rather one 
flux quantum. In the original basis, the vector potential 
A falls off as l/r for distances greater than a penetra- 
tion depth away from a single vortex. The transformed 
vector potential. A, follows the superfluid velocity Vs(r), 
and vanishes exponentially for distances greater than a 
penetration depth. 

The matrix model of the previous sections, Hj, uses 
as a basis the eigenstates of the full Hamiltonian TYq, 
which includes the magnetic field. The models described 
by Hjj and Hjjj first apply the above transformation, 
then use the the basis defined by the eigenstates of just 
the T-invariant part of the Hamiltonian. This has two 
important effects when the magnetic field is weak. 

First, weak T-breaking implies that the phase gradi- 
ent of the order parameter varies slowly on the scale of 
a coherence length: £,qV9 ^ 2tt. Therefore the term in 
the bare Hamiltonian due to the superfiuid velocity is 
small compared to the energy gap: (l/m)p • eA(r) ^ 
Vp\79 <C Aq. Second, weak T-breaking in a conven- 
tional superconductor means that the deviation in the 
magnitude of the order parameter from the bare BCS 
value is small: Sf/Ao = /(dr/17)|/(r) - AqI/Aq < 1. 
In the Ginzburg-Landau description of an isolated vor- 
tex, the order parameter relaxes to the bare value ex- 
ponentially, on the scale of the coherence length. Hence 
Sf/Ao ^ Co/^H — H/Hc2- In a more realistic descrip- 
tion, /(r) may have power-law behavior |/(r) — Ao| ~ 
(^o/r)" which yields Sf/Ao - {H/Hc2f/^ [n = 1), or 
Sf/Ao - {H/Hc2)HHa2/H) {n = 2). In either case, 
rotating to a basis of the eigenstates of the T-invariant 
part of Tio for which the order parameter is real yields 
a Hamiltonain which is a small correction to the BCS 
theory. 

The models Hjj and Hjjj are defined by the assump- 
tion that the matrix elements of the T-breaking correc- 
tions to the BCS theory are uncorrelated. That is, in 
Eqs. (|39|) and @, the matrices M and M^^) (the super- 
fluid velocity correction), and A/^^^ (the order parameter 
correction), are assumed to be random. The difference 
between the two models is that in Hjj the effect of the 
deviation of the order parameter from the constant mag- 
nitude is ignored, whereas in Hjjj it is included. Hence in 
Hji the pairing is diagonal, AqSij, because {i\T\j) — Sij 



(T is the time reversal operator T : ipi{r) ip*{r)), 
while in Hjjj, the deviation of the order parameter yields 
an additional correction: {i\Tf{r)\j) — AoSij + a2Ml-' , 
where 02 ^ Sf. 

These models describe small corrections to BCS and 
hence properly reproduce the BCS result as H/Hc2 0. 
The matrix model of the previous sections, Hj, can not 
be described as a small correction to BCS. In the weak 
field hmit, W/Ao ~ H/Hc2 ^ 0, so /i(w) ^ 5(w). How- 
ever, because the bare level spacing is always assumed to 
be much less than the T-breaking scale {Sq ^ W), the 
resulting matrix model still pairs many levels with equal, 
random amplitude. This leads to large deviations from 
BCS. For W/Aq 0, the density of states has a node 
near zero energy, p{E) ^ E, rather than a gap as in the 
BCS spectrum Q. 

At weak fields, or more generally for weak T-breaking, 
when the correction to BCS behavior can be consid- 
ered small, the most appropriate random matrix descrip- 
tion is that of the third model listed above: H 



in- 



ks 

the T-breaking energy scale increases, this description 
must break down. Although the order parameter can 
always be written as the BCS value plus a correction, 
/(r) = Aq -I- Sf{r), the assumption that the matrix ele- 
ment {i\T Sf{r)\j) is random is no longer justifiable when 
Sf/Ao ^ 1; /(r) and Sf{r) are comparable functions. At 
sufficiently high fields, the only consistent assumption of 
randomness applies to the entire pairing matrix, as as- 
sumed in Hj. This may called 'chaotic' pairing. 

Figure 1 illustrates schematically the different regimes. 

Physically, at low fields, vortices are well separated, 
bound states are confined primarily to individual vor- 
tices, and the effect of the superfiuid velocity on extended 
states is small, yielding only a small correction to the or- 
dered BCS state. At higher fields, core states overlap 
significantly and form energy bands, which interact in 
a complicated way with extended states, the superfluid 
velocity, and whatever disorder is present. In the Ander- 
son metal-insulator transition, delocalization is accom- 
panied by a change from ordered (Poisson) statistics to 
chaotic (random matrix) statistics. In a type-II super- 
conductor, we are suggesting that the delocalization due 
to increasing vortex density (increasing magnetic field) 
is accompanied by a transition from a small correction 
to the BCS theory to chaotic pairing. This raises inter- 
esting possibilities and questions, discussed in the next 
section, regarding the nature of such a transition in a 
superconductor. 



VI. SUMMARY AND CONCLUSIONS 

When the disorder and magnetic field are sufficiently 
large in a type-II superconductor, so that electronic 
structure calculations and small deviations from the BCS 
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theory are no longer applicable, the chaotic pairing pic- 
ture of the random matrix model may apply. The moti- 
vation for the model is that the density of states is not 
sensitive to details of the rapidly varying 0(1) factors in 
the matrix elements between electrons in a Cooper pair, 
but rather to the average structure of the pairing matrix. 

In section III we showed that the average amplitude to 
pair electrons is related to the Cooperon response func- 
tion previously discussed in the context of weak localiza- 
tion theory. For a dirty metal this response function may 
be calculated and yields a simple form for the average 
pairing amplitude, characteristic of a diffusive system. 

The connection between the average pairing amplitude 
and diffusion is related to what was called by de Gennes 
"ergodic pairing". That connection was established by 
de Gennes in the context of evaluating the critical field 
of dirty superconductors, Hc2, using perturbation theory. 
Here we have shown that this connection exists indepen- 
dent of perturbation theory. This allows us to evaluate 
the density of states for fields far away from the regime 
in which perturbation theory is valid. 

Finally, we considered the weak field limit of the ran- 
dom matrix model, H ^ Hc2- At sufficiently weak fields 
the random matrix picture must break down because the 
system must approach the ordered BCS state. This sug- 
gests a crossover between a low-field ordered state, in 
which the corrections to BCS theory are small and can 
be understood as being due to states localized near in- 
dividual vortices, to a high- field chaotic state, in which 
most of the states are extended. A useful analogy is the 
Anderson metal-insulator transition. Instead of varying 
the electron density (chemical potential) from a regime 
of localized states to a regime of delocalized states, in 
the superconductor we vary the vortex density (magnetic 
field) between a regime of bound states to a regime of ex- 
tended states. In the metal-insulator transition we have 
a change in statistics from Poisson to random matrices; 
in the superconductor we have a change in the pairing, 
from BCS-like to chaotic. 

The connection between the field-dependence of the 
pairing in a superconductor and the nature of the low- 
energy quasiparticles raises several interesting questions. 
Is there a true localization phase transition for quasi- 
particles? Is there a sharp change between an ordered 
regime and a chaotic regime or merely a crossover? Are 
there simple models for which, using techniques similar 
to those in the theory of disordered metals, the chaotic 
pairing picture may be derived? What is the interplay be- 
tween localization, disorder, and magnetic field for high- 
Tc superconductors, where the prediction in zero field 
for a two-dimensional d-wave superconductor is "univer- 
sal conductivity" [Q? A further important issue is to 
test the validity of the random matrix prediction for the 
density of states, for each of the possible regimes, by di- 
rect comparison with tunneling experiments. 
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FIG. 1. In the absence of impurities and for an ordered 
lattice of vortices, the quasiparticle spectrum of a type-II su- 
perconductor may be calculated by numerically solving the 
Bogoliubov-de Gennes equations. In the presence of random 
disorder, and for weak magnetic fields, the quasiparticle spec- 
trum may be calculated as a small correction to the BCS 
theory, Eq. (^^, similar to the Abrikosov-Gor'kov theory. 
For stronger disorder or magnetic fields, the system enters 
a chaotic regime in which the pairing is completely random, 
Eq. (H). 
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